Multiple solutions of CCD equations for PPP model of benzene 
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To gain some insight into the structure and physical significance of the multiple solutions to the 
coupled-cluster doubles (CCD) equations corresponding to the Pariser-Parr-Pople (PPP) model of 
cyclic polyenes, complete solutions to the CCD equations for the A^^ states of benzene are obtained 
by means of the homotopy method. By varying the value of the resonance integral /3 from —5.0 eV to 
—0.5 eV, we cover the so-called weakly, moderately, and strongly correlated regimes of the model. For 
each value of /3 230 CCD solutions are obtained. It turned out, however, that only for a few solutions 
a correspondence with some physical states can be established. It has also been demonstrated that, 
unlike for the standard methods of solving CCD equations, some of the multiple solutions to the 
CCD equations can be attained by means of the iterative process based on Pulay's direct inversion 
in the iterative subspace (DIIS) approach. 
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I. INTRODUCTION 

The understanding of the electronic correlation effects 
in cyclic polyenes (or [Mjannulenes, with the chemical 
formula CmHm, where M — AniQ + 2, mo = 1,2,...) has 
attracted considerable interest for many years. To some 
extent this interest has been caused by the fact that these 
molecules provide models for extended one-dimensional 
metallic-like systems which epitomize the difficulties en- 
countered in the description of extended many-electron 
systems. A very useful tool in these, studies proved to 
be the coupled cluster (CC) method.El Even its simplest 
variant, the CC method with double excitations (CCD)jll 
is able to reproduce most of the electronic-correlation ef- 
fects. More sophisticated CCSD0 and CCSD(T)kl meth- 
ods provide often quite accurate approximations to the 
full configuration-interaction (FCI) results (which define 
the limit for a given orbital basis set). The CC method 
furnishes a set of coupled nonlinear inhomogeneous equa- 
tions, of unknowns representing the amplitudes of the CC 
operator T (the t amplitudes), in which the number of 
equations is equal to the number of unknowns. However, 
due to their nonlinearity, the CC equations have multi- 
ple solutions; thisijDroblem was first analyzed by Zivkovic 
and Monkhorst.oB It is still little known about the math- 
ematical prQfierties of these solutions. 

RecentlyjElO for the first time, some insight into the 
structure of the complete set of solutions to the CC 
equations has been gained by aonlying the powerful 
homotopy (continuation) methodll^'U to the CCD and 
the CCSD equations corresponding to some four-electroH. 
four-orbital systems, known as the H4 and P4 models.t3 
These studies included the correspondence Jpetween CC 
and configuration interaction (CI) methods,!] and the in- 
fluence of the approximate form of the cluster operatat 
on the structure of solutions of the associated equations.^ 
Within the symmetry and spin-adapted CCD formalism. 



the CCD equations for the H4 model comprise a system 
of six coupled quadratic inhomogeneous equations for six 
unknowns (the t2 amplitudes of the double-excitation 
operator T2). According to the Bezout theorem, the 
maximum number of solutions amounts in this case to 
2^ = 64. Yet Kowalski and JankowskiQ showed that the 
CCD equations for the H4 model case have only 12 solu- 
tions (6 real, and 3 pairs of complex solutions). Some 
of these solutions corresponded to the variational re- 
sults of the configuration-interaction with doubles (CID) 
method. 

In the present paper we study, by means of the homo- 
topy method, multiple solutions of the CCD equations for 
the six-electron six-orbital-system corresponding to the 
Pariser-Parr-Pople (PPP)E3T3 model of benzene. This 
system is of special interest for theoretical chemists: the 
sextet of TT-electrons is responsible for the peculiar chem- 
ical properties of benzene which is the archetypal aro- 
matic molecule. The PPP model is known to provide|-a 
simplified yet reliable picture of many-electron systems.t£l 
For benzene, one gets an almost perfect description of the 
TT-electron part of the excitation spectrum by optimizing 
the ampirical parameters of the PPP model at the FCI 
level JiZlli3 Moreover, in the PPP model the strength of 
the electronic-correlation effects can be controlled with- 
out changing the molecular geometry, by simply adjust- 
ing the value of the so-called resonance integral f3 (< 0). 
Very high symmetry of the PPP benzene model ensures 
that the number of symmetry-independent parameters 
entering the exact (FCI) or an approximate (e.g., CCD) 
wave function is small. All that makes the PPP benzene 
model particularly suitable for testing quantum-chemical 
methods dealing with the electronic correlations. 

Benzene is the first member of the [Mjannulenc fam- 
ily for which, as mentioned above, the understanding of 
the electronic structure is of rather general significance. 
These systems, especially when described within the PPP 
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model, seem deceptively simple: their one-electron states 
are fully determined by the spatial symmetry, and so 
is their Hartree-Fock (HF) determinantal wave function. 
Yet [MJannulenes turn out to be a challenge to the exist- 
ing computational methods of treating the electronic cor- 
relations: it has been found that for larger [A/]annulcnes 
(corresponding to mo > 2), in the so-called strong corre- 
lation regime of /3, the CCD method breaks down com- 
pletely, as no real solution of the CCtLeauations exists 
for P greater than some critical valueolHj A recent CC 
study of [MJannulenes, taking into account the double 
(D), triple (T), and quadruple (Q) excitations in the CC 
operator, showed that even the CCDTQ method breaks 
down for these systems jssihen the correlation effects be- 
come sufficiently strong.Eil For the PPP benzene model, 
a real CCD solution representing the ground state can 
be found for any /3 < 0. However, as f3 approaches 0, 
the CCD t2 amplitudes deviate more and more from the 
corresponding FCI values, thus suggesting that some cor- 
relation effects peculiar to the higher [MJannulenes may 
surface already in the strongly correlated regime of the 
PPP benzene model. This provides an additional moti- 
vation for the present study. 



The point-symmetry group of benzene is Dgh, but its 
subgroup Ce is sufficient for the symmetry considerations 
in the 7r-electron approximation. The molecular orbitals 
of the TT symmetry (ttMOs), expressed as linear combina- 
tions of ttOAOs, are completely determined by the pro- 
jections onto the irreducible representations of the Cg 
group and read as 



(1) 



where k — 0,±1,±2,3 stand for the symmetry labels. 
For k 0,3 the ttMOs are complex, ip^ = V'-fe, and, due 
to the time-reversal symmetry, correspond to degenerate 
orbital energies, e(fc) — e(— fc). In the restricted Hartree- 
Fock (RHF) description of the ground state of benzene 
the occupied ttMOs correspond to fc = 0, ±1 and the 
unoccupied ttMOs to fc = ±2, 3. The formulas for the 
orbital energies and the HF 7r-electron energy read asEJ 

£{k) = 2/3cos(fc7r/3) + 7° - [v{k) + v{k + 1) + v{k - 1)J , 

(2) 



II. PPP MODEL OF BENZENE 

A detailed description of the PPP model of annulcnes 
may be found in Rcf. In benzene the C atoms form a 
regular hexagon, and the C-C bonds are assumed to be 
of the length = 1.4 A. The PPP model invokes the tt- 
electron approximation and describes the six tt electrons 
of benzene by using a minimal basis set of 2pz atomic 
orbitals associated with the six carbon atoms. The set 
of these tt atomic orbitals (ttAOs) is then subject to tka 
symmetrical orthonormalization procedure of L6wdin,E2l 
yielding the set of six orthonormalized tt atomic orbitals 
(ttOAOs), denoted by Xm , m = 0, ±1, ±2, 3 (we use here 
the numbering system employed in Ref. 21). The Fock- 
space Hamiltonian H for benzene, built according to the 
prescriptions of the PPP model, is given in Eq. (2) of 
Rcf. |2l|. The following semiempirical parameters are used 
in the PPP model: 

(i) a, the so-called Coulomb integral for the carbon atom, 
representing the binding energy of electron described by 
the ttOAO Xm (for simplicity we put a = eV), 

(ii) /3(< 0), the resonance integral, corresponding to the 
electron transfer between the neighboring ttOAOs, Xm 
and Xm+i (our results are for /3 = —5.0, —4.0, —3.0, 
-2.5, -2.0, -1.5, -1.0, and -0.5 eV), 

(iii) 7mn, representing the two-center two-electron in- 
tegrals (xmXnIXmXn); usually One calculates = 
l{Rmn), where Rmn is the distance between the cen- 
ters of orbitals Xm and Xn and function 7(i?) is given 
by some simple analjytical formula. We use the Mataga- 
Nishimoto formulae^, 7(-R) = e^[R + e'^{'-i°)'^]^^ , where 
e is the electron charge and 7° — 7(0) = 10.84 cV. 



Shf = 8/3 + 37° - [3w(0) + Av{l) + 2v{2)] . 



(3) 



These formulas depend on parameters /3, 7*^, and the 
two-electron integrals calculated in the ttMO basis: 

^ 6-^ [7" + 2cos(g7r/3)7(i?°) 

+ 2cos{2qTr/3)j{V3R°) + cos(g7r)7(2i?")J , (4) 

where q — 0,1, 2, 3, and the modulo-6 rule of addition is 
assumed for the symmetry labels. 



III. CCD METHOD FOR THE PPP MODEL OF 
BENZENE 

In the standard single-reference CC theory, the 
ground-state FCI wave function ^E* for a six-electron sys- 
tem is represented as 



^' = exp(r)$. 



(5) 



where T = Ti+T2-\ h Tg is the CC operator and $ is 

the RHF determinantal wave function, playing the role 
of the reference configuration. The T„ components of the 
CC operator correspond to the connected n-tuple exci- 
tations from occupied to unoccupied spin-orbitals. Each 
Tn operator depends on some linear parameters, here- 
after referred to as the tn amplitudes; the ordered set 
of all tn amplitudes form a vector denoted by t„. In 
the CC theory one introduces a similarity-transformed 

Hamiltonian, H = exp(— T)_ff exp(T), whose amplitudes 
hPf - are certain connected functions of the amplitudes 
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of the Fock-space Hamiltonian, and the t„ amphtudes. 
The electronic correlation energy for the ground state ^' 
may be calculated in the CC theory as 



^corr = - ^HF = /l(tl,t2) - ^HF 



(6) 



i.e., it is a function of only ti and t2 amplitudes. The t„ 
amplitudes {n — 1, 2, . . . , 6) may be calculated by solving 
the set of CC equations: 



/lf(tl,t2,t3) 
^i'i'(tl,t2,t3,t4) 



= 0, (7a) 
= 0, (7b) 



TABLE I. Correspondence between indices of Eq. (^) and 
quasimomentum indices fci, k2, and q. 

i ki k2 q 



($ 



abcdef 
ijklrnn 



abcdef 
ijklmn 



(ti,t2,t3,t4,t5,t6) = 0, (7c) 



where 



where ^fj*- - is an n-tuply excited configuration. Writ- 
ten in an explicit form, the CC equations (|^) form a set 
of coupled inhomogeneous nonlinear equations, with the 
number of unknowns (the t„ amplitudes) equal to the 
number of equations. On a basis of the FCI method, it 
can be shown that the exact tn amplitudes are real. In 
the simplest approximate variant of the CC method, the 
CCD one, one puts T — T2 and neglects ii, ^3, and ^4 
amplitudes in Eqs. (^) and (7b) (as well as the remaining 
CC equations). The CCD equations ( [7b| ) then become 
a set of coupled inhomogeneous quadratic equations for 
the unknown ^2 amplitudes, and from Eq. (^ an approx- 
imate correlation energy, E^^^ , is calculated. 

In the PPP model of benzene the occupied and un- 
occupied orbitals belong to different representations of 
the Ce group, which causes Ti and T5 to vanish by sym- 
metry. Because ti =0, the RHP function $_becomes 
equal to the Brueckner determinantal function.cJ In this 
case Eqs. ( fz^ ) are automatically satisfied, and the CCD 
method becomes equivalent to the CCSD one. When the 
nonorthogonal spin adaptation of the CCD equations is 
performed (see, e.g., Ref. ^), the spin-adapted t2 ampli- 
tudes for benzene may be written as i(fci, k2, q), where ki 
and ^2 are occupied-TrMOs labels, and q (> 0) is chosen 
such that ki + q and k2 + q are unoccupied- ttMOs labels 
(see Ref. ||). It can be shown that there are 11 differ- 
ent sets of fci,A:2, and q, thus defining 11 ^2 amplitudes 
for benzene. By assuming that these amplitudes are real 
and employing the tiarie-reversal symmetry, one finds a 
symmetry constrainiEU 



t{ki,k2,q) = t{-k2,~ki,q) 



(8) 



which reduces the number of the symmetry-independent 
t2 amplitudes for benzene to 8. 

In Ref. ^ we employed the general non-orthogonally 
spin-adapted CCD equations for the PPP model of the 
[MJannulenes, with the t2 amplitudes <(fci, fc2, q) subject 
to the symmetry constraint (^. For benzene, the set of 
these equations may be written as 



8 8 



(9) 



1,2,. 



and the unknowns 



J 



1,2,..., 8 stand for the symmetry-independent ^2 am- 
plitudes, see Table |. The formula for the electronic- 
correlation energy now reads as 



E. 



CCD 



(10) 



The linear coefficients in Eqs. (|^) and (pO|): at, bij (for 
* 7^ j): Cijfc = Cikj, and dj can be expressed as some 
linear combinations of the two-electron integrals defined 
in Eq. (Q). Only the diagonal elements ba depend on 
parameter /3: 



Ae, - 2v{0) , 



(11) 



where 

Aci = Ae(A;i, /C2, q) = e{ki + q) + e{k2 - q) - e(fci) - e(fc2) , 

(12) 

through the dependence on f3 of the orbital energies (||) . 
There is still some symmetry hidden in the set of CCD 
equations (|^): due to the so-called alternancy symmetry 
of the PPP Hamiltonian H (sec Ref. and references 
therein) one finds that 



X2 



t(0, 0, 2) = t(l, -1,2)^X3^ 



(13) 



(for the general formula for [MJannulenes, see Ref. |T]). 
This property propagates into Eq. (^) making them in- 
variant with respect to the interchange of indices 2 and 
3. However, our set of CCD equations is not explic- 
itly adapted to the alternancy symmetry, and solutions 
breaking this symmetry are, in principle, possible. 

The usual method of solving the CCD equations is 
based on an iterative procedure, which in the case of 
Eq. (O) may be written as follows: 



j=i k=i 



= -{Aay^ [ai + ^(&y - AeiSij)x 

8 8 

+ ^ ^ CijkXj Xf. ] , 

j=i fc=i 



(n) 
j 



(14) 
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with xf^ = 0, i = 1,2,.. .,8. The first iteration fur- 
nishes t2 ampUtudes that substituted into Eq. (|l^) give 
the second-order M0Uer-Plesset (MP2) correlation en- 
ergy. When convergent, this simple iterative procedure 
provides a single solution to the CCD equations (^. Such 
a solution is bound to be real and to preserve the alter- 
nancy symmetry, which corresponds to the fulfillment of 
Eq. (|l|). In our calculations for the PPP model of ben- 
zene with P e [—5.0 eV,0 eV], we found that the above 
described simple iterative procedure is indeed conver- 
gent, and the correlation energy calculated by substitut- 
ing the convergent t2 amplitudes into Eq. (|l^) approxi- 
mates the FCI value for the ground state. The agreement 
between the CCD and FCI results {t2 amplitudes and 
Ecorx) is very good in the weakly and moderately corre- 
lated regimes (/? in the vicinity of —5.0 eV and —2.5 eV, 
respectively), but becomes rather poor in the strongly 
correlated regime (/3 > —0.5 eV). The CCD and FCI re- 
sults for (i = —2.5 eV and /3 = —0.5 eV may be found in 
Ref. M. 



IV. MULTIPLE SOLUTIONS OF CCD 
EQUATIONS FOR BENZENE 

The CCD equations (||) comprise a set of 8 cou- 
pled quadratic inhomogeneous equations (with real co- 
efficients) for 8 unknowns. According to the Bezout the- 
orem, such equations may have up to 2^ — 256 solutions, 
complex in general. In principle, a complete set of solu- 
tions can be fwp4 by means of the homotopy (continua- 
tion) method.Hllil Below we present the results obtained 
by applying the homotopy method to the CCD equations 
(^ corresponding to various values of the resonance inte- 
gral /3: from —5.0 eV (representing the weakly correlated 
regime) to —0.5 eV (representing the strongly correlated 
regime). The FCI resjttlts used for comparison were cal- 
culated with GAMESS.E3 

Equations (^ have been derived by taking into account 
the spin and the time reversal symmetries, as well as the 
reality of the t2 amplitudes, see Eq. (||). However, these 
equations may have also complex solutions. Such solu- 
tions have to appear in pairs: if x = (xi, X2t ■ . ■ ,xg,) is a, 
complex solution, then its complex-conjugate x* is also a 
solution. For a complex solution, the correlation energy 
calculated from Eq. (|l^) assumes (in general) a complex 
value, and the complex-conjugate solutions correspond to 
the complex-conjugate values of the correlation energy. 
Pairs of complex solutions, x and x*, will be called de- 
generate, since they correspond to the same real part of 
the complex correlation energy calculated from Eq. (10). 
Some solutions may also violate the equality ( |l3|) derived 
from the alternancy-symmetry. Such symmetry-broken 
solutions, for which X2 7^ 2:3, also have to appear in 
pairs: if x is such a solution, then x', in which the val- 
ues X2 and are interchanged, has to be a solution as 
well. Pairs of the real symmetry-broken solutions, x and 



x', are also degenerate. In the case of the general com- 
plex symmetry-broken solutions, x and x*, x' and x'* 
form a degenerate quadruplet. We have found, however, 
a special class of complex symmetry-broken solutions in 
which the only complex values correspond to X2 = x'^. 
In such a case one has x* = x', and the pair of solu- 
tions x and x' corresponds to the same real value of the 
correlation energy calculated from Eq. (|o|). Thus, the 
solutions of Eqs. (^) may be classified into five distinct 
categories: real symmetric (which are non-degenerate), 
real symmetry-broken, complex symmetric, general com- 
plex symmetry-broken, and a special class of complex 
symmetry-broken with the real energy. 

In Table || the number of CCD solutions belonging 
to different categories is presented for several values of 
the resonance integral /3. The total number of solutions 
(230) is surprisingly large: it is only slightly smaller than 
the Bezout upper bound (256), and much larger than 
the number of solutions for the 114 model (6 real and 

6 complex, compared to the upper limit 64 allowed by 
the Bezout theoremH) . The number of CCD solutions is 
also much larger than the number of the FCI solutions 
for the PPP model of benzene, having the same symme- 
try as the RHF wavefunction $, see further discussion. 
This implies that most of the CCD solutions have no 
physical meaning. Since the FCI method is equivalent 
to the full coupled-cluster (FCC) method, it is interest- 
ing that the CCD equations, which have smaller order 
and a smaller number of the unknowns, have more solu- 
tions than the more sophisticated FCC equations. This 
suggests that the truncated CC equations are unable to 
utilize all the symmetries that are present in the many^ 
electron Hamiltonian. For the 114 model it was observedil 
that the number of solutions increased from CCD (12 so- 
lutions) to CCSD (60 solutions), and then decreased to 

7 for FCC. 

While for the 114 model the number of pairs of complex 
solutions was equal to the number of real solutions, for 
the PPP model of benzene there are much more complex 
than real solutions. Since all the FCI results are real, 
the complexity of the solutions must be caused by the 
truncation of the CC operator. Since all the linear coef- 
ficients in the CCD equations (^) are real, each complex 
solution must have its complex conjugate counterpart. 
However, while solving these equations by using the ho- 
motopy method, we have found in several instances that 
some complex-conjugate solutions were missing (dupli- 
cated solutions were also encountered). Since the prob- 
lem has not been previously reported, we think that it 
may have arisen due to the numerical complexity of the 
problem (a lot of roots to be traced in the homotopy al- 
gorithm) . It brings in the question whether the CCD so- 
lutions presented here are complete. After removing du- 
plicates and adding solutions that must be present due to 
symmetry, we have found that the total number of states 
is constant (equal to 230) for each considered value of (3. 
The number of solutions preserving the alternancy sym- 
metry has also appeared to be constant (equal to 124). 
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TABLE II. Numbers of states of different symmetry for various (3 given in eV. 
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It seems therefore unlikely that certain states has been 
overlooked. 

Due to the alternancy symmetry of the PPP model, 
the six-electron symmetry states of benzene ^ht into 
two categories, denoted by "minus" and "plus."Ea Includ- 
ing the spin symmetry and the spatial symmetry of the 
Dgh group, the ground-state RHF wave function <i> cor- 
responds to the symmetry label Aig. Among 400 FCI 
states of benzene generated by GAMESsEII there are only 
18 states corresponding to the A'^^ symmetry. These 
TO^Jg-states (m = 1,2,..., 18) are, in general, non- 
orthogonal to and thus may be expressed in the form 
of the CC expansion (||). For each case, the correspond- 
ing tn amplitudes may be extracted from the FCI linear 
coefficients: there are 7 ^2 amplitudes, 2 amplitudes, 
7 amplitudes, and only a single amplitude (only 
non-redundant parameters are counted). 

We are interested in identifying solutions to the CCD 
equations which have physical significance, i.e., which 
correspond to some states of the ^A^^ characteristic rep- 
resented in the model by relevant solutions to the FCI 
equations. In order to gauge a similarity between the 
^2^*-''^ amplitudes corresponding to a given CCD solution 
and the ^2*^^ amplitudes corresponding to some FCI so- 
lution of the A^g symmetry, we use parameters 9 and 7] 
defined below: 



kCCD+FCI 



tCCD||tFCI| 



(15) 



where we use the real part of the amplitudes t2^^ ■ Here 6 
measures the angle between the vectors, and rj — the ratio 
of the vector lengths; in the above analysis we use the vec- 
tors corresponding to the full set of t2 amplitudes (of the 
dimension 11), i.e., containing the symmetry-redundant 
amplitudes fulfilling Eq. (H). In principle, a complex so- 
lution may be considered an approximation to some real 
solution of the FCI equations as long as the imaginary 
parts of the ^2 amplitudes are small in comparison to the 
real parts. 

Some of the CCD solutions, obtained for several val- 



ues of /3, are characterized in Table III. Solution numbers 
JT-CCD are assigned in accordance with the increase (of the 
real part) of the corresponding correlation energy value 
calculated from Eq. (|l^). For instance, tt-ccd = 005 de- 
notes the 5th solution. Due to the large number of solu- 
tions, we consider only those that are the closest to some 



FCI ones, i.e. those corresponding to 9 and ry closest to 
and 1, respectively. (The complete set of solutions may 
be obtained from the authors). 

Solution 005 for each of the /3 values (except for (3 = 
—4.0 eV, where it-ip, solution 007) is considered to be 
the ground state,llaEl and can be obtained by applying 
the standard iterative process of Eq. (|lj). Indeed, this 
solution is the most similar to the FCI ground state Aig 
for all the (3 values studied (in the strongly— correlated 
region the similarity is, however, rather pooiEj). Except 
for that state, there is little similarity between the CCD 
solutions considered and the FCI states, both in energy 
and amplitudes. 

For all values of (i the state corresponding to the 
ground state has, in fact, the energy larger than some 
of the other CCD solutions. We have found a pair of real 
symmetry-broken "underground" solutions that have the 
correlation energy ranging from —265.333498 eV for (3 = 
-5.0 eV to -165.644985 eV for (3 = -0.5 eV. Moreover, 
for each (3 there is a pair of complex symmetric solutions 
(except for (3 — —4.0 eV, where there are two pairs of 
such solutions) with the real part of E^^^ lower than the 
ground-state CCD value. For instance, for (3 = —0.5 eV 
one has E'^^^ = -13.907202 ± il2. 763757 eV, and for 
(3 = -5.0 eV E'^^^ = -3.892588 ± illO.045436 eV. All 
the "underground" solutions do not resemble any of the 
FCI states; such solutions have not been observed for H4 
model. 

Since the standard method of obtaining CCD solutions 
is via iterative procedure (14), it may be useful to study 
the the performance of this procedure in attaining multi- 
ple solutions when starting from various CCD solutions 
obtained by the homotopy method. Obviously, if the so- 
lutions were exact, the iterations would stop after the 
first iteration. However, the small errors in the numeri- 
cal values would lead to a non-trivial iteration sequence. 
The question is whether this sequence is convergent and, 
if yes, what is the converged result. We tested all the 
real solutions and it is surprising that the series is either 
divergent or the result is identical to the ground state. 
It means that only the ground state is stable in the iter- 
ative process. It explains why the iterative process with 
the MP2 starting point gives only the ground state solu- 
tion. 

A different behavior has been found when applying 
the dirpc.t inversion in the iterative subspace (DIIS) 
methodnSEj. It turned out that this method, which is 
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TABLE III. Correspondence between some states of the PPP model of benzene and solutions to the CCD equations for 
various /3 values. FCl and CCD correlation energies are in eV, nccD is CCD solution number, (in rad) and rj parameters are 
defined in Eq. @. 
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useful for accelerating convergence, improves also the sta- 
bility of the iterative process. Most of the states are sta- 
ble in the DIIS iterative process and only between 6-8 
states out of the 22-32 real ones either diverge or con- 
verge to a different state (not necessarily to the ground 
state). It is worth stressing that the "underground" so- 
lutions are iteratively unstable for all the cases tested. 

V. CONCLUDING REMARKS 

In this paper we have investigated the complete set of 
solutions of the CCD equations corresponding to the PPP 
model of benzene. This is the largest system for which 
the complete set of solutions has been obtained. We have 
found that the number of solutions (for a broad range of 
the resonance integral /3) is surprisingly large (equal to 
230), approaching the limit given by the Bezout theorem 



(equal to 256). One may wonder if some peculiar prop- 
erties of the PPP model are responsible for that prolif- 
eration of solutions. To this end, we have checked the 
behavior of the PPP model of butadiene, which is anal- 
ogous to the lid, model studied previously by Jankowski 
and Kowalski.lZna We have found, however, that the CCD 
equations for these two models have roughly the same 
number of solutions. We thus conclude that the large 
number of solutions of the CCD equations correspond- 
ing to the PPP model of benzene have to be related to 
some special properties of a cyclic 7r-electron system. The 
number of solutions is expected to grow exponentially for 
the larger cyclic polyenes (annulenes). The known diffi- 
culties with solving the CC equations for these systems 
(see Ref. ^ and references therein) are undoubtedly a 
consequence of this multiple-solution problem. On the 
other hand, some non-standard solutions of the CCD 
equations for larger annulenes may turn out to be similar 
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to certain excited-state FCI solutions of these systems. 
Since attaining these solutions by means of the homotopy 
method seems to be out of reach in the near-iuture, we 
find it encouraging that the DIIS methodcJ'E^ of carry- 
ing the iterations in the CC method proved quite effective 
in assessing the stability of several multiple solutions to 
the CCD equations for benzene. This has inspired us to 
look for multiple solutions of the CCD equations for the 
next member of the annulene family, cyclodecapentaene 
(CioHio). Indeed, several new solutions were found by 
combining a sort of random generation of the initial ^2 
amplitudes with the DIIS iterations. The results of these 
investigations will soon be published. 
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